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Abstract. A new non-disruptive error field (EF) assessment technique not restricted 
to low density and thus low beta was demonstrated at the EXTRAP-T2R reversed field 
pinch. Stable and marginally stable external kink modes of toroidal mode number 
n=10 and n=8, respectively, were generated, and their rotation sustained, by means 
of rotating magnetic perturbations of the same n. Due to finite EFs, and in spite 
of the applied perturbations rotating uniformly and having constant amplitude, the 
kink modes were observed to rotate non-uniformly and be modulated in amplitude. 
This behavior was used to precisely infer the amplitude and approximately estimate 
the toroidal phase of the EF. A subsequent scan permitted to optimize the toroidal 
phase. The technique was tested against deliberately applied as well as intrinsic error 
fields of 7i=8 and 10. Corrections equal and opposite to the estimated error fields 
were applied. The efficacy of the error compensation was indicated by the increased 
discharge duration and more uniform mode rotation in response to a uniformly rotating 
perturbation. The results are in good agreement with theory, and the extension to 
lower 77, to tearing modes and to tokamaks, including ITER, is discussed. 
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1. Introduction 

Error Fields (EFs) [lj are known to lower confinement in toroidal plasmas and, in 
general, to lower their rotation. An exception is the Neoclassical Toroidal Viscosity 
induced "offset rotation" [2]. Additionally, EFs can seed magnetic islands by EF 
penetration or cause locking of pre-existing islands [31 H] , often resulting in disruptions 
[5]. Correction techniques were developed to minimize EFs and their effects pp. The 
most widespread technique consists in reducing the density until a locked mode forms, 
with the achievement of lower densities indicating better EF correction: performing the 
ramp in presence of different, deliberately applied static magnetic perturbations (MPs), 
allows to indirectly infer the EF and its best correction. This, however, is the best 
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correction at low density and thus at low beta, but experimental evidence [6] and recent 
ITER modeling [7] suggest that the best correction at high beta can be significantly 
different due to plasma response [7j. Furthermore, the technique described above needs 
several discharges, each terminating with a locked mode and, often, a disruption. 

Here we present an EF correction (EFC) technique not based on the appearance 
of low-density Locked Modes, not restricted to low beta, non destructive, and requiring 
only a fraction of a discharge. In fact, it could be deployed multiple times within the 
same discharge, to dynamically assess and correct the EF in real time as it evolves. 

The technique was inspired by DIII-D experiments in which rotating MPs were 
used to rotate and reposition locked tearing modes of poloidal/toroidal mode number 
m/n=2/l and bring them in view of gyrotron launchers, for their stabilization [8]. In 
spite of the MPs rotating uniformly, the experiments showed non-uniformities in the 
mode rotation. These observations were explained with the fact that rotation was 
taking place in the presence of EFs, and the mode was in fact locked to the resultant 
of the static EF and the applied rotating MP. Such resultant does not rotate uniformly, 
neither is it constant in amplitude. In fact, if the MP is weaker than the EF, the 
resultant doesn't even describe complete revolutions and, for short periods of time, it 
can even rotate opposite to the MP [8]. 

The idea behind the present work is that (1) the applied MP is known, (2) the 
EF+MP resultant can be measured from the dynamics of a mode (not necessarily a 
tearing mode) locked to it and therefore (3) the EF can be deduced from the difference. 
Specifically, we applied uniformly rotating MPs and analyzed the non-uniform mode 
rotation and the variations of mode amplitude to precisely infer the amplitude and 
approximately estimate the toroidal phase of the EF. A subsequent scan permitted 
to optimize the toroidal phase. The experiments were carried out at the EXTRAP- 
T2R reversed field pinch [U] and were made possible by the coil control system briefly 
described in Sec. 2. For simplicity, and unlike the earlier DIII-D experiments with 
m=2, n=l locked tearing modes, the modes used here were external kink modes of 
m=l and high toroidal numbers, n=8 and 10. The rationale for this decision, explained 
in Sec. 3, is that at EXTRAP-T2R these modes are stable, which made the experiment 
easier. Sec. 4 contains numerical predictions of the time evolution of the amplitude and 
toroidal phase of a stable kink mode rotating in the presence of a resistive wall, static 
EFs and a rotating MP. The predicitions are based on a time-dependent cylindrical 
MHD model for thin-wall boundary conditions [10] and are in good qualitative and 
quantitative agreement with the experimental results presented in Sec. 5. The transition 
from incomplete to complete mode rotation allowed to quantify the strength of the n=8 
and 10 EFs. Corrections of equal strength and varying toroidal phase were applied by 
means of static MPs until the toroidal phase was optimized. In this optimization we 
used two indicators of good EFC: the extended discharge duration and a more uniform 
mode rotation in response to a uniformly rotating perturbation. 
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2. Experimental Set-up 

EXTRAP-T2R |9| is a reversed field pinch (RFP) of major radius i?=1.24m and plasma 
minor radius a=0.183m. External to the vessel is a double- layer Copper shell of thickness 
d—2x 0.5mm. A magnetic field of poloidal mode number m and toroidal mode number 
n penetrates across this resistive shell (or wall) in a time T mn which is a fraction of 
the "long" resistive wall time t w . Theoretically, in the cylindrical limit, this evaluates 
t w = fJ,ocrr w d— 13.6ms, where a and r w are the conductivity and radius of the wall. 
Experimentally, r w was measured to evaluate approximately r w =l 1.2ms, slightly larger 
than earlier estimates [11]. 

Earlier works [T2l [T3] provided evidence of EFs at various n. EFs at n=+8 and 
n—+10, which the present paper focuses on, might be due to the vacuum vessel being 
made of bellows sections joined with short flat sections (15 in total, placed at irregular 
toroidal angles). 

EXTRAP-T2R is equipped with a system for the active control of resistive wall 
modes (RWMs) and EFs. External to the vessel, and just inside the shell, is an array 
of 4 (poloidal) x 32 (toroidal) saddle loops used to sense the radial field, depicted in 
Figjlj Located outside the shell is an array of 4x32 control coils or actuators, powered 
by suitably modified audio amplifiers, used to apply radial field perturbations. Various 
control algorithms were developed and tested with success [HI [15], but of particular 
relevance here is a feedback algorithm presented in [IB] , capable of simultaneously 
sustaining preset finite amplitudes for all modes of m=l and —16 < n < 15. These 
amplitudes are typically set to 0, to suppress RWMs and EFs. 

The algorithm was modified for the sake of the experiments presented here: the 
capability was added to suppress all modes except 1-4 modes of up to 4 different n. Those 
unsuppressed modes were not controlled in feedback. For the corresponding values of 
n, it was possible to apply static MPs of given amplitude and toroidal phase and/or 
rotating MPs of given amplitude, frequency and phase. In particular, rotating MPs were 
used to excite stable external kink modes and drive their rotation. On occasions, static 
MP were used to deliberately apply "proxy" EFs, stronger than the intrinsic machine 
EFs, thus easier to detect and optimal to test the method validity. 

3. Discussion on Modes and Mode Numbers 

The RFP equilibrium used for this experiment is defined by the reversal parameter [T7] 
F = B^(a) I {B^ = —0.25 and the pinch parameter = B g (a) / (5^)=1.66, where Bg 
is the poloidal field and (B^) is the toroidal field B^ averaged over the plasma cross- 
section. Fig. [2^, summarizes the stability (growth rate 7 <0) and instability (7 >0) 
of m—1 modes of various toroidal numbers n. The diagram was computed for the 
equilibrium used. However, small changes of equilibrium can have significant effects, 
thus the value of the diagram is mostly illustrative. Modes of m 7^1 are typically stable 
for the RFP equilibrium used (low G) and will not be considered here. The convention 
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on the sign of n is as follows: modes with a handedness corresponding to the pitch of 
field lines inside the reversal surface have n < and modes outside the reversal surface 
have n >0. The reversal surface is defined as the surface where changes sign, located 
at approximately r/a=0.82. 

Stable modes can only exist in presence of finite perturbations of the same m and 
n. They form (disappear) on a timescale 7 _1 as soon as those perturbations are applied 
(removed) . 

Note that "modes" here can refer to tearing or kink modes, depending whether the 
corresponding rational surface (resonance) is internal to the plasma and an island can 
form, or not. If external to the plasma, the perturbation can only kink the plasma itself. 
Very high positive n, though, can resonate at a surface external to the reversal surface, 
but inside the plasma. 

From left to right in Figj2]we have: 

(i) rapidly rotating (tens of kHz) Tearing Modes (TMs) of n <12. These modes are 
internal, resonant (q = —m/n in the plasma) and stable |48] . 

(ii) unstable and marginally unstable kink modes (—11 < n < 6, except n = —1, which 
calculations suggest being marginally stable) growing on a time scale of the order 
of r w or slower. These are basically RWMs [9]. 

(iii) stable kink modes (n > 7) only exisiting in presence of an MP of the same n. These 
mode grow (decay) on timescales of order t w or slower as soon as the MP is turned 
on (off). However, they do not grow indefinitely in time: they only grow up to a 
saturated value proportional to the MP strength. The proportionality factor, or 
gain, is dictated by the plasma response. These modes can also be referred to as 
stable RWMs [13]. 

Correspondingly, the following considerations can be made with regard to suitability 
to EF assessment by driven rotation: 

(i) Shielding by the shell makes it impossible to couple rotating MPs to fast-rotating 
modes like the TMs mentioned above. Besides, previous work showed that, at 
EXTRAP-T2R, the main effect of static MPs and EFs on TMs is their magnetic 
braking [18] . MPs rotating at few tens of Hz are expected to have the same effect. 

(ii) Unstable and marginally unstable RWMs can in principle be used, but their growth 
would be in competition with other effects on their amplitude, making it difficult 
to extract the effect of the EF. These modes will be the subject of future work. 

(iii) RWMs of n > 8, however, are stable. This simplifies the extraction of EF and 
MP effects on the RWM amplitude, compared with unstable RWMs of lower n 
mentioned above. Also, being stable, they decay as soon as their MP+EF drive 
is zeroed, which is also an indication of good EFC More generally, as soon as the 
drive is reduced or steered, so is the RWM. Finally, stable RWMs amplify EFs, 
making it easier to measure them with sensor coils. 
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It should also be mentioned that a suitable mode for the technique proposed does 
not need to be preexistent: it can be generated when necessary, for instance by a /3 
ramp at low rotation, like NTMs in the DIII-D experiments cited above |B]. It does not 
need to be rotating either. For example, it can be a non-rotating locked Tearing Mode 
not preceded by a rotating precursor, but rather seeded by EF or MP penetration. On 
the other hand, it does not need to be exactly locked either: a slowly rotating Quasi- 
stationay Mode will still be affected by EFs in its rotation. 

Rather, the mode of choice needs to interact with the EF. The EF and, as a result, 
the mode, needs to depend on the toroidal angle 0, and it must be possible to define 
a ^-dependent potential energy, reaching a minimum for a certain value of <fi. For this, 
it is not necessary for the mode and the dominant EF to have the same m and n: the 
interaction of modes and EFs of different mode numbers will still depend on cf). 



4. Numerical Predictions 



Let us expand the radial field b r measured at the wall in poloidal and toroidal 
components b™ ,n in the cylindrical (large aspect ratio) limit: 



b r (r, 6, 0, t) = b ? ,n ( r > *) exp(im9 + 



(1) 



where r, 9 and <f> are the radial, poloidal and toroidal coordinate respectively. 
In the absence of perturbations and control the growth rate r y m ' n of a RWM of mode 
numbers m and n is defined by 

d 



dt 



b m,n _ m.njm.n = q 



(2) 



In presence of a thin shell (thin compared with the skin depth for the mode and 
timescale of interest), 7 m,n is related to r w ~1 1.2ms by [3] 



7 T w 



d{rb™> n ) 
Or 



(3) 



The quantity in square brackets represents the jump of the derivative across the thin 
shell and r~ and r+ denote radial locations just inside and outside the radial location 
of the wall (and, with good approximation, of the sensors), r w . 

In presence of a time-dependent external perturbation b™^, Eq. [2j describing the 
time-evolution of the radial-field harmonics evaluated at the wall, 6™' n , modifies as 

follows [31 HOI H9l ETil EH E21 E31 H21 El]: 



d 



,m,n lt«,» 
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'ext 



dt r r 
where we adopted the notations of 
for b™^™, related to r w by 



~d(rb v r ) 
dr 



(4) 

and T m,n is the wall diffusion time 



(5) 



Error Field Assesment from Driven Mode Rotation 



6 



Here b v r is the radial component of the vacuum field, i.e. in absence of plasma, but 
in presence of the thin resistive wall. The jump in the radial derivative across the shell 
on the right hand side of the equation is due to eddy currents induced in the conductive 
thin shell itself. These eddy currents can be due to modes in the plasma and/or currents 
in the control coils. 

The n dependence of r y m ' n and T m,n in EXTRAP T2R for m=l for the plasma 
equilibrium used in the present work are plotted in Figj2j 

In the remainder, for simplicity, superscripts m,n will be omitted from &™' n , 6^.'", 
r m ' n and 7 m ' n . For our purposes b ext consists of a static field of amplitude b (the 
resultant of the intrinsic EF, proxy EF, if any, and deliberately applied static MP, if 
any) and a rotating field of amplitude bi\ 

K^^bo + he^, (6) 

where u is the angular frequency at which the field, or equivalently the currents in 
the active coils, oscillate. A mode of toroidal number n rotates with angular velocity 
u/n. 

Note that both 6 an d b\ are complex and determine the phase of b ext (t) and that 
Eq. [i] can be analytically solved. Assuming as initial condition b r (0) = 0, we obtain: 

&r = — (e* - 1) + -. - , (e* - e*"*) (7) 

This expression describes the behavior of the radial magnetic field perturbation 
measured at the sensors. This measured field (Eq|7j) is due to the applied external 
perturbation (Eq|6]) but is modified by the plasma response. 

An initial transient of order 7 _1 is visible in Figj3j This transient typically lasts 
few ms (see plots of 7 in Figj2^i, for various modes). After the transient, Eqj7]for stable 
modes (7 <0) simplifies as follows: 



b r = \b r \e iwt+A * (8) 



where 



\K\ = -7== 2 0) 

A(p = arctan — (10) 

7 

Eqj7|is used to simulate the experimental results with r w = 11.2ms and r ) m ' n and 
T m ' n from Fig|2j Fig|3] shows simulation results for n=10, a static field 6 = 0-4mT and 
four values of the probing field b\, rotating at 50Hz. Note that in the experiment b\ is 
known, while b is the unknown total static EF, which can be determined as follows. 

A small rotating field 61 b^r^Juj 2 + 7 2 does not drive a complete rotation of the 
mode, but only an oscillation around the phase of b (here assumed 0), accompanied by 
large amplitude variations (blue curves in Figpl). Under the effect of a large rotating 



Error Field Assesment from Driven Mode Rotation 



7 



field bi 3> b T^oo 2 + 7 2 , instead, the sensors register a complete, uniform rotation and 
smaller relative changes of amplitude (green curves). Finally, when the rotationally 
shielded probing-field magnitude — , bl (Eq.9) is slightly smaller or larger than the 

Ty/u> 2 +'y 2 

EF b , the result is a large oscillation in the first case (red), and a complete, although 
highly non- uniform rotation in the second (black). Correspondingly the angular velocity 
has negative or positive peaks, respectively (Fig(3fc). The transition from one behavior 



(oscillations) to the other (complete rotations) occurs for — , bl = bo. Experimentally 

identifying this transition allows to indirectly measure the EF amplitude bo- Also, 60 
has toroidal phase 0=0 in the example of Fig|3j but this can be easily generalized: if its 
phase is 0, the measured field will oscillate around 0. can also be measured when b r 
reaches its minimum, i.e. as the mid-point of the large phase oscillation recorded before 
the mode describes a complete rotation. However, that is a large oscillation occurring 
in a short time, hence this measurement of suffers from errors. In fact, the closer one 
gets to the transition, the preciser is the measurement of |6 |, but the higher is \d(f>/dt\, 
thus the more rapid the variation of and the less precise its measurement, as if the 
two quantities |&o| and were related by an indetermination principle. 

Fig |3] displays curves for four values of the rotating field b±, but the calculations were 
repeated for several other values. The results are summarized in Figjl], as a function 
of b\. In particular, Fig|4^ shows the minimum the minimum sensor measurement 
that can be reached at some time (i.e., for some toroidal orientation of the applied MP) 
for a given amplitude b\ of the rotating probing field. Fig|4]D presents the modeled phase 
velocity of the measured perturbation, d<p 1,10 /dt, reached at that same time. As already 
seen in Fig|3j when bl' 10 reaches a minimum, the angular velocity reaches an extreme 
(either a negative minimum, if the rotation is incomplete, or a positive maximum, if the 
rotation is complete). 

Perfect error field cancellation is obtained when the modeled measured b], ,l ° is 
suppressed to zero (b± ^0.6mT in Figj4^i). At that time, the angular speed d<f) l ' 10 /dt 
transitions from strongly negative to strongly positive (Fig|4](b)). 

Note that observing a transition for a certain value of b\ in Figsj3]j4] does not imply 
that the EF 60 takes that same value. This would only be true in the limit of very slow 
rotation, uj — )■ 0. Otherwise, the amplitude b\ of the rotating perturbation (as applied 
outside the shell) undergoes some screening, and a stronger b\ needs to be applied 
to cancel a given bo- From Eqj9} we can estimate the screening factor to evaluate 
\/t\Juj 2 + 7 2 ~ 0.7 for w/27r=50Hz. Thus, observing a transition at b\ ~0.6mT in 
Figs(3]j4] agrees with the fact that the EF assumed above for this example evaluates 
&o=0.4mT. 



5. Experimental Results 



Here we present the proof of principle of the proposed EF assessment technique based 
on driven rotation of modes. For the reasons discussed in Sec. 3, it was decided to use 
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stable external kinks and characterize EFs of the same mode numbers: m=l, n=8 and 
m — 1, n=10. We started with n=10 kinks (more stable) and progressed to n=8, which 
was found to be the lowest stable n and to have a smaller, thus, more difficult detect, 
intrinsic EF associated with it. 

5.1. Results for n= 10 

For the first test, static MPs were used to apply "proxy" EFs of n=10, i.e. fictitious 
EFs, stronger than usual and thus easier to detect. The n=10 EF was suppressed by 
automatic feedback as usual [16l dH [15] in the first 10ms of the discharges illustrated 
in Fig|5j At t >10ms, though, the intrinsic n=10 EF was left uncontrolled and an MP 
was applied, to simulate an n=10 EF of 0.92mT. In order for such EF to be constant, 
the actual programmed MP applied was increasing slightly on a timescale comparable 
with the discharge duration. This timescale is unrelated to 7~ 1 , t w or r m,n , but 
rather due to subtleties in the control hardware 

In the same time- interval (t >10ms) an n=10 rotating MP was also applied. Its 
rotation frequency, 50Hz, enables the observation of at least two rotation periods within 
a typical EXTRAP-T2R discharge. The resultant of the static MP, rotating MPs, and of 
the much smaller uncorrected EF, excites a RWM of n=10 and thus of negative growth 
rate (Figj2^i), that is, stable. The time-evolution of its amplitude, toroidal phase and 
angular velocity is regulated by EqjTjor, after a transient of few ms, by its approximation, 
Eqj8] Indeed we obtained sensor coil measurements (Figsj5|6]) very similar to the 
calculations in Figsj3|4j including the hyperbolic behavior of d^ 1 ' 10 / dt. In particular, the 
transition from incomplete to complete mode rotation occurs between &!=0.35mT and 
6i=0.55mT (Fig|5j3) , at b\ = 0.48 ± 0.05mT (Figj^). For consistency and comparison 
with EqjTj we retained in the sensor measurements the direct contribution of the applied 
MPs (both static and rotating). In other types of analysis, e.g. to extract the field b r 
associated with a mode, the contribution of the MPs could have been removed by 
"vacuum subtraction" , i.e. by taking separate measurements in presence of MPs but in 
absence of plasma, and subtracting them from the actual plasma measurements. This 
was not the case here, because b and b\ defined above include the applied static and 
rotating MPs. On the other hand, the effect of the vertical (VF) and toroidal field 
(TF) coils on the sensor measurements was vacuum-subtracted. The reason is that 
those coils introduce a high background in the b r measurements, but such background 
is nearly perfectly axisymmetric, i.e. the VF and TF are not significant sources of n=10 
EF. 

With the screening factor 1/Va/w 2 + 7 2 taken into account, it is "estimated" 
that the proxy EF (which is actually known, but is treated as an unknown for this 
test) amounts to bo = 0.89 ± O.lOmT. The error results from the propagation of the 
uncertainties in the measured screening factor (circa ±10%) and in the rotating MP 
amplitude (while the programmed amplitude was constant, the actual amplitude was 
observed to fluctuate by approximately ±5%). The bo = 0.89 ± O.lOmT measurement 
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agrees with the applied field evaluating 6 =0.92mT although, rigorously, it needs to be 
compared with the EF+MP resultant, which is done below. 

The toroidal phase can be estimated from Figj5^i to be ~ (0 ± 0.1)7r, as this is 
the value around which signals oscillate. Again, this is in agreeement with the applied 
proxy having 0=0. Alternatively, the EF phase could be inferred from the measurement 
of at the time when b r reaches its minimum i.e. when the EF (the proxy EF) is best 
corrected by the EFC (the rotating MP). However, care should be exerted because 
varies rapidly in the vicinity of that optimal value (except of course if the rotating MP 
is strong enough to dominate the kink dynamics and make it rotate uniformly). 

The screening factor is l/r^/u 2 + 7 2 , and in Secji it was calculated to evaluate 
approximately 0.7, on the basis of the prescribed to, calculated 7 (Figj2^i) and calculated 
r (Fig^p). In the experiment, however, rather than calculating it, we estimated 
the screening factor experimentally by comparing sensor measurements in two plasma 
discharges with, respectively, large static and large rotating MPs. The MPs were large 
in order for the intrinsic EF to be negligible. Comparing the two discharges showed 
that, for m—1 and n=10, rotating MPs of a;/27r=50Hz generate, in presence of plasma, 
weaker signals (by a factor 0.54±0.05) than static MPs of the same intensity. 

For a precise determination of 0, as well as to confirm the good quality of the EFC, 
corrections were applied, of amplitude 0.89mT and various 0, in two sets of discharges, 
and two indicators of good EFC were monitored. The first indicator is the duration of the 
discharge: it is well-known that an improved EFC leads to better confinement and thus 
longer discharges, as indeed shown in Figj7^i. The broad maximum at = (1.05±0.25)7r 
agrees with the expected best EFC having = ir. 

In another set of discharges, the EFC phase was scanned from one discharge 
to the other, a small, uniformly rotating MP was applied, and the uniformity of the 
measured mode rotation was studied: for perfect EFC, the mode was expected to rotate 
uniformly as well. Fig|7)D confirms that the smaller variation of angular velocity during 
a rotation period is obtained for = (0.9 ± 0.1)7r, in agreement with expectations. 

To further illustrate the good EFC, Figj8] shows a reconstruction of the plasma edge 
and radial field based on actual sensor data. Two features can be clearly recognized: 
the strong m—1, n=10 perturbation at a time when the rotating MP has the same 
as the proxy EF, and its minimization half a rotation-period later (when the MP and 
proxy EF nearly cancel each other). 

After additional tests with a smaller proxy EF (0.23mT, not shown for brevity), 
proxy EFs were eventually removed and the the technique was applied to the assessment 
of the intrinsic m—1, n—10 EF. The results are presented in Fig|9]and analyzed in FigfTOj 
The best correction is obtained for a rotating perturbation &i=0.11mT, corresponding 
to a static (unshielded) b = 0.20±0.02mT, in agreement with earlier estimates [T2"| [T3]. 
Figj9]D indicates that its toroidal phase is approximately 0=0. 

In retrospect the total static field error EF+MP in the proxy experiment presented 
above amounted to 1.12 ± 0.04mT and had phase ~ (0 ± 0.1)7r, agreeing within two 
standard deviations with the measurements, 60 = 0.89 ± O.lOmT and ~ (0 ± 0.1)7r. 



Error Field Assesment from Driven Mode Rotation 



10 



5.2. Results for n=8 

The n=10 results presented above are easily extended to n >10 at EXTRAP-T2R, as 
the corresponding RWMs are even more stable. For this reason, it is more interesting 
to extend the technique to lower n. The extension to unstable n will be the subject of a 
future work. In the present Section, we concentrate on marginally stable modes. It was 



experimentally determined that the lowest positive stable n is n=8 (Fig. 11), in rough 
agreement with the prediction in Figj2^i, that modes of n > 7 are stable. 

Similar to the n=10 experiments, a probing n=8 MP rotating at 50Hz was applied, 
and its amplitude was varied from one discharge to the other. In a first set of discharges 



(Figs 12 13), an n—S, 0=0 proxy EF of 0.41 mT was applied. The screening factor 
was 0.57 ± 0.06 and the transition to complete rotations occurred at &i=0.44mT. This 
implies a bo = 0.77 ± 0.08mT estimate which should not be compared with the static 
MP, but with the EF+MP resultant (see below). The toroidal phase is <f)/ir = 0.1 ± 0.1. 
Ultimately the technique was applied in the absence of proxy EFs, to assess the 



intrinsic m=l, n=8 EF (Figs, 14 15). It was established that it has a toroidal phase 
<f)/n = 0.5 ±0.1 and an amplitude of 0.15±0.02 mT, in agreement with earlier estimates 
[j~2| IT5] . A scan of the discharge duration confirms that the best EFC is obtained 



for a phase (j}/ir = 1.3 ± 0.2 (Fig. 16), which is opposite to the EF phase 4>/ir = 0.5, as 
expected. 

With this information taken into account, it is deduced that the vector sum of the 
intrinsic EF and applied static MP in the n=8 proxy experiment described above had 
amplitude bo = 0.44mT, agreeing within two standard deviations with the measurement 
bo = 0.77 ± 0.08mT, and orientation 0/-7T = 0.11 ± 0.1, in good agreement with the 
measurement, <f)/ir = 0.1 ± 0.1. 

Summary and Discussion 

In summary a new non-disruptive error field (EF) assessment technique not restricted 
to low density and thus low beta was demonstrated at the EXTRAP-T2R reversed field 
pinch. Stable Resistive Wall Modes (RWMs) of toroidal mode number n=8 and 10 
were generated and their rotation sustained by rotating magnetic perturbations. Due 
to finite EFs, and in spite of the applied perturbations rotating uniformly and having 
constant amplitude, the RWMs were observed to rotate non-uniformly and be modulated 
in amplitude (FigsjSJ [6j |9j [lOj [l2j 13 14 15,). This behavior was interpreted with a 



simple theoretical model (Eq(7|) and used to characterize and correct intrinsic n=8 and 
10 EFs, leading to longer discharges and more uniform mode rotation (Figsj7, [16). 

The probing perturbation was rotated (i.e., its phase was scanned) during the 
discharge, while its amplitude was scanned from shot to shot. It is conceivable that 
in another device with sufficiently long discharges, such as ITER, both parameters 
(amplitude and phase) can be scanned within a single discharge or even a fraction of it, 
e.g. by deploying a rotating, growing ( "spiraling" ) probing field. While fast relative to 
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the discharge duration (to guarantee several rotation periods and a thorough scan), for 
ease of interpretation the MP rotation and growth should be slow compared with the 
RWM growth, the resistive wall time and current diffusion time. 

The technique works best with non-growing modes driven by or sensitive to EFs 
(or, more precisely, to the EF+MP resultant). The principle is that when the drive is 
reduced or rotated, the mode shrinks or rotates accordingly. At EXTRAP-T2R we used 
stable Resistive Wall Modes (RWMs) of n=8 and 10 for the reason that these modes, 
indeed, do not grow. There was no fundamental reason behind the high n. 

The technique is easily extended to classical or neoclassical tearing modes in other 
devices, e.g. n=l locked modes at DIII-D, due to EF penetration or not, provided they 
naturally saturate, as it is often the case for periods of several hundreds of ms, or are 
kept small by stabilizing Electron Cyclotron (EC) Heating. EC Current Drive is even 
more stabilizing but its effect depends on its toroidal phase relative to the island O-point 
and, in fact, can become destabilizing if deposited in the X-point. This can be used to 
the advantage of the technique presented, as it enhances the variations of the measured 
field, but attention has to be paid to separate the effects of the EF and of the current 
drive. 

The extension to unstable modes such as n = — 1 RWMs at EXTRAP-T2R is left 
as future work. 
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Figure 1. Schematic of (from inner to outer) vacuum vessel, sensor coils, resistive 
shell and control coils at EXTRAP-T2R. 
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Figure 2. (a) m=l RWM growth rate and (b) wall diffusion time as a function of n 
for the equilibrium used in the present work. 
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Figure 3. Calculated (a) amplitude, (b) phase and (c) angular velocity of m=l, n=10 
component of radial field at the sensor coils, modeled after EqjTJ for EF amplitude 
b = OAmT and MP rotating frequency w/27r=50Hz. The amplitudes of the rotating 
perturbations are b\ = 0.2mT (blue curves), bi — 0.6mT (red), b\ = 0.9mT (black) 
and b\ = 2mT (green). Note that the static EF evaluates 0.4mT, but the transition 
occurs at a rotating field amplitude 0.6mT< b\ <0.9mT. 
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Figure 4. (a) minimum amplitude of 6^ 10 as calculated in Fig{3^i, as a function of 
the MP strength 6 X for EF amplitude b — 0.4mT and MP rotation frequency 50Hz. 
(b) corresponding phase velocity at the time of minimum amplitude. Error bars are 
smaller than symbol sizes. 
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Figure 5. Experimental (a) amplitude, (b) phase and (c) lms-smoothed angular 
velocity of m=l, n=10 component of radial field at the sensor coils, in presence of 
"proxy" EF of amplitude bo = OAmT and phase 4>=0 (same as n Figj3|. MPs rotate at 
w/27r=50Hz, and have various amplitudes (see Figj6]). Transition to complete rotation 
occurs between b\ =0.35 mT and bi =0.55 mT. Signals are only plotted for the duration 
of the discharges. Error bars are comparable with or smaller than symbol sizes. 
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Figure 6. (a) minimum amplitude of &*' 10 evaluated from Fig[5^ and alike, as a 
function of the MP strength bi for EF amplitude 6o=0.92 mT and MP rotation 
frequency 50Hz (same values as in Fig|4|. (b) corresponding phase velocity at the 
time of minimum amplitude. 




Figure 7. Optimization of the toroidal phase <fi of the n=10 EFC and evidence of 
good EFC for optimal <fi, as evinced from (a) longer duration of the discharge and (b) 
more uniform rotation. 




Figure 8. Reconstruction of m=l, n—10 external kink mode (stable RWM) based 
on actual saddle coil measurements at times when the MP (a) reinforces or (b) nearly 
cancels the EF. For clarity the radial deformation is exaggerated by a factor 10. 
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Figure 9. Like Fig[5j except that no "proxy" EF is applied. Here the only static bo is 
the intrinsic m=l, ri=10 EF. Transition to complete rotation occurs between bi =0.078 
mT and bi =0.116 mT. 
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Figure 10. Like Fig|6j but for the n=10 intrinsic EF results of Fig(9] 
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Figure 11. Measured amplitude of m—1 modes of various n, in discharges where the 
corresponding n was left uncontrolled at i=30-60ms. The n—7 signal, unstable, is due 
to a growing RWM destabilized by the uncontrolled n—7 EF. Other signals (n=8, 9 and 
10) take finite but small, non-growing (stable) values, due to intrinsic EFs, possibly 
amplified or attenuated by the plasma response. Note that both the EFs and plasma 
response can fluctuate with time. 




Figure 12. Similar to FigjSJ but for n=8. Also, feedback is disabled and the static and 
rotating perturbations applied at t >20ms instead of t >10ms. Transition to complete 
rotation occurs between b\ =0.42mT and b\ =0.46mT. 
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Figure 13. Like FigM but for the n—8 results of Fig 12 




Figure 14. Like Fig[l2| except that no "proxy" EF is applied. The only static bo is the 
intrinsic m=l, n=8 EF. Transition to complete rotation occurs between b\ =0.08mT 
and bi =0.095mT. 
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Figure 15. Like Figfl3| but for the intrinsic n—8 results of Fig 14 




Figure 16. Optimization of toroidal phase </> of intrinsic ?i=8 EF correction (EFC). 
Longer discharge durations are indicative of better EFC. 



